A Complete Basis for a Perturbation Expansion of the General 

A^-Body Problem 

W. Blake Laing*, David W. Kelle^ Martin Dunn, and Deborah K. Watson 
Homer L. Dodge Department of Physics and Astronomy, University of Oklahoma 

(Dated: February 19, 2009) 

Abstract 

We discuss a basis set developed to calculate perturbation coefficients in an expansion of the 
general A^-body problem. This basis has two advantages. First, the basis is complete order-by-order 
for the perturbation series. Second, the number of independent basis tensors spanning the space for 
a given order does not scale with N, the number of particles, despite the generality of the problem. 
At first order, the number of basis tensors is 23 for all N although the problem at first order scales 
as A^^. The perturbation series is expanded in inverse powers of the spatial dimension. This results 
in a maximally symmetric configuration at lowest order which has a point group isomorphic with 
the symmetric group, Sm- The resulting perturbation series is order-by-order invariant under the 
A^! operations of the Sm point group which is responsible for the slower than exponential growth 
of the basis. In this paper, we perform the first test of this formalism including the completeness 
of the basis through first order by comparing to an exactly solvable fully-interacting problem of N 
particles with a two-body harmonic interaction potential. 
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I. INTRODUCTION 



In a previous paper[l], we described the development of a perturbation method for the 
general S-wecve A'"-body problem through first order. Group theoretic and graphical tech- 
niques were used to describe the interacting A^-body wave function for a system of identical 
bosons with general interactions. Solutions for this problem are known to scale exponentially 
with N requiring that resources be essentially doubled for each particle added [2, 3]. As A'" 
increases beyond a few tens of particles, this growth makes a direct numerical simulation 
intractable without approximations given current numerical resources. Typical approxima- 
tions truncate the Hilbert space of the exact solution by using a basis that spans a smaller 
Hilbert space or by truncating a perturbation series or both[4-14]. With bosonic systems, 
various Monte Carlo approaches may be employed which scale in a polynomial fashion with 
N , making larger- calculations feasible. [7-9, 11] 

In Paper 1, a perturbation series is developed in inverse powers of the spatial dimension. 
This results in a maximally symmetric configuration at lowest order having a point group 
isomorphic to the symmetric group Sn . The basis used is complete, at each order finite, 
and, in fact small, having only twenty three members at first order, despite the scaling 
of the problem at this order. 

In this paper, we test this formalism which truncates the perturbation series, but deter- 
mines each term in the series exactly using group theory and graphical techniques. This 
perturbation series is order-by-order invariant under the operations of the ^'jv point group. 
The tensor blocks needed at each order can thus be decomposed into a basis also invariant 
under the Sn point group. It is this restriction, i.e. the invariance under N\ permutations, 
that slows the growth of the number of basis tensors as increases, resulting in a basis that 
is small in contrast to the A^^ growth of the vector space at first order. 

We have named the basis tensors, "binary invariants." "Binary" because the elements 
within a basis tensor are ones or zeros; "invariants" because the basis tensors are invariant 
under N\ permutations of the particles labels. 

This strategy effectively separates the A^ scaling away from the rest of the physics and 
then tackles the A^-scaling problem using the symmetry of the Sn group. The full problem, of 
course, scales exponentially in A^, so as higher orders of the perturbation series are included 
the full exponential A^ scaling of the problem will appear. However, in this methodology, the 
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N scaling problem has been compartmentalized away from the rest of the physics and dealt 
with using group theoretic and graphical techniques; i.e. it becomes a straight mathematical 
problem. Once this mathematical work, which involves significant analytical effort, has been 
completed at a given order, it never has to be repeated again for a new interaction or a 
different number of particles, 

The formalism being tested in this paper has been presented in a series of papers beginning 
with the isotropic, lowest-order ground-state wavefunction [15], the isotropic, lowest-order 
ground-state density profile [16], and the isotropic, first-order ground-state wavefunction [1]. 

In Section II we review the large-dimension point group symmetry and discuss its impli- 
cations in perturbation theory. In Section III we discuss the binary invariants themselves 
and the graphs which label them, and in Section IV we very briefly review the general theory 
for the wave function through first order from Paper 1. In Section V we compare the results 
of the general theory for the wave function through first order applied to the problem of the 
harmonically interacting system under harmonic confinement with an expansion of the exact 
interacting wave function through first order. This work is performed in Appendix A where 
we exactly solve for the wave function of the harmonically interacting A'^-particle problem 
under harmonic confinement in D dimensions. This exact solution is then expanded through 
first order in the square root of the inverse dimension of space. 

II. THE LARGE-DIMENSION, POINT GROUP SYMMETRY AND IT'S IMPLI- 
CATIONS IN PERTURBATION THEORY 

As discussed in previous papers on dimensional perturbation theory (DPT) the Hamil- 
tonian and Jacobian-weighted wave function and energy are expanded in powers of 6^^"^ : 
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where 



(2) 
(3) 
(4) 



The internal displacement coordinate vector yl is a column vector composed of the displace- 
ment radii and the displacement angle cosines 7^^: 
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and is related to the dimensionally-scaled internal coordinates by 

n^f^ + S'/^f', and 7y =7oo + <^'^'%-, 



(7) 



where the fj are the dimensionally scaled radii and the jij are the angle cosines between the 
position vectors of the particles. 

The superprescript on the F and G tensors in parentheses in Eqs. (3), (4) and (5) denotes 
the order in S^^"^ that the term enters (harmonic being zeroth order). The subprescripts 
denote the rank of the tensors. 

In general, iJ„ is of order n + 2 in the elements and derivatives of y' (2"^^ order in the 
derivatives), and formed of cither all even or all odd powers of the elements and derivatives 
of y' when n is even or odd respectively. 
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According to Eq. (7) the system localizes as D — > oo on a configuration centered about 
= ^oo and = 7^. This structure has the highest degree of symmetry where all 
particles are equidistant from the center of the trap and equiangular from each other (a 
configuration that's only possible in higher dimensions). The point group of this structure is 
isomorphic to Sn which in effect interchanges the particles in the large dimension structure. 
This together with the fact that the full D-dimensional Hamiltonian, H , is invariant under 
particle exchange means the DPT expansion of Eq. (1) is order-by-order invariant under 
this Sn point group, i.e. the Hj are each invariant under the point group. This greatly 
restricts the F and G tensors of Eqs. (4) and (5). In three dimensions a corresponding 
A^-particlc structure would have a point group of lower symmetry, i.e. one not isomorphic 
to Sn despite the fact that all of the particles are identical. It is this profound restriction 
on the F and G tensors from the point group symmetry, which in itself is a direct 
consequence of developing a perturbation theory about the large-dimension limit, that allows 
for an essentially analytic solution of the problem at a given order in the perturbation 
theory for any N . For example, with a third-rank tensor, instead of (A'"(A'" -|- l)/2)^ possible 
independent elements (a number which becomes quickly intractable since, for example, with 
only ten particles this would be more that a million elements), the S]\f point group symmetry 
restricts the number of independent elements to only twenty three, independent of the value 
of A^. 

A. The Reducibility of the F and G Tensors under Sn 

The maximally symmetric point group , together with the invariance of the full Hamil- 
tonian under particle interchange, requires that the F and G tensors be invariant under the 
interchange of particle labels (the Sn group). In fact the various blocks of the F and G 
tensors are themselves invariant under particle interchange induced by the point group. For 
example, ''3 -P'/^JIj.j) is never transformed into ^^3 F^i]) (^ki){mn) ' ^ labels are 

preserved since the group does not transform an coordinate into a 7'^^^-) coordinate. 

The various r — 7 blocks of the F and G tensors may be decomposed into invariant, and. 
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this time irreducible, blocks. Thus for example, ^"^Q may be decomposed into the blocks 

i < j , k < I 

all of which remain disjoint from one another under the group. Significantly, invariance 
under particle interchange requires that tensor elements related by a label permutation 
induced by the point group must be equal. This requirement partitions the set of tensor 
elements for each block into disjoint subsets of identical elements. Consider the elements of 
the ^^2Q^^ block. The element ^2'*Qi[i belongs to a set of N elements (of the form ^2'*Qi,i) 
which are related by a permutation induced by the point group, and therefore must have 
equal values. Likewise, the element ^^2QY2 belongs to a set of A?" — 1 elements related by 
a permutation induced by the point group and sharing a common value. Proceeding in 
this fashion, we observe that the blocks of the lowest-order tensors are partitioned into the 
following set of identical elements which remain disjoint under particle interchange: 

V iandA; (9) 

V i^j8.ndk^l (10) 
ijim) = ^?QSpn) V i<k, Z < m , and n > p (11) 
IJmn) ^ j <k andly^rrKn (12) 
JkiUki) V i < J and A; < Z (13) 
Vm),{in) ^ i<j^k>i and l<m^n>l (14) 
JmnUpq) ^ i^jj^kj^landm^u^p^q (15) 

i < j , k < I m < n , p < q 

The lowest-order block matrices contain the sets of elements in Eqs. (9)- (15) arranged in an 
intricate pattern. A similar, but more involved partitioning occurs for higher rank F and 
G tensors. 
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In Paper 1 it was shown that this decomposition could be expressed in terms of binary 
tensors which are invariant under Sn and are labelled by graphs. We term these binary 
tensors, binary invariants. 

III. BINARY INVARIANTS 

A. Introducing Graphs 

Definition 1 A graph Q = (V, E) is a set of vertices V and edges E. Each edge has one 
or two associated vertices, which are called its endpoints/i?/. 

For example, . . is a graph Q with three vertices and three edges. We allow our 
graphs to include loops and multiple edges [18]. A graph contains information regarding the 
connectivity of edges and vertices only: the orientation of edges and vertices is insignificant. 

We introduce a mapping which associates each tensor element with a graph as follows: 

1. draw a labeled vertex (• i) for each distinct index in the set of indices of the element 

2. draw an edge ( LJ) for each double index {ij) 

3. draw a "loop" edge (O « ) for each distinct single index i 

For example, the graph corresponding to the tensor element ^^2Q7lij) under this mapping is 

Two graphs with the same number of vertices and edges that are connected in the same 
way are called isomorphic. The elements of the group arc permutation operations which 
interchange particle labels. Two elements with graphs that are not isomorphic are never 
related by a permutation of Sn group. We label each set of isomorphic graphs by a graph 
with no vertex labels. Denoting the set of unlabeled graphs for each block as GxiX2...x„, 
where n is the rank of the tensor block (and therefore the number of edges in each graph in 
the set) and X is r or 7 , we have 

Grr = {00, 0} 

G^r = {0-,0 -} (16) 
= {O, J,Z} 
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Gr = {0} 

G, = M 

Grrr = {9^,8^^0 0} (17) 
G,rr = {g., 00,0-0, 00,00} 

Gj^r = {C»,00,>^,I-0'LO,rP,~0} 

<^jjj = {e,^,o-,_L>U>2>I],s} 

Now consider a tensor, for which all of the elements labeled by a a single isomorphic set of 
graphs are equal to unity, while all of the other elements labeled by graphs heteromorphic to 
this single set of isomorphic graphs are equal to zero. We term this tensor a binary invariant, 
[B{Q)]^^^,^2,--- since it is invariant under the Sn group, and we label it by the graph Q; sans 
particle labels at the vertices, for the non-zero elements all of which are equal to unity. Thus 
each of the above graphs denotes a binary invariant, [-B(^)]i/i,;/2,... ■ Explicit expressions for 
these binary invariants for arbitrary N may be found in the EPAPS document [19]. 

B. Binary invariants are a basis 

That the small number of binary invariants are a complete basis with which to represent 
any invariant tensor in the N{N + l)/2-dimensional f'-7' space can be seen as follows. 

Lemma 1 The set of binary invariants, [-B(^)],/i,i^2,- > f^''^ ^^^^ ^ ^'^^ linearly independent 

Proof This follows from the fact that no binary invariant shares a non-zero element with 
another binary invariant for a different graph. 

Lemma 2 The set of binary invariants for all Q & G spans the invariant tensor space. 

Proof From the mapping of Sec. Ill A above, relating graphs to tensor elements in the 
r'-7' space, every tensor element is related to an unlabeled graph, and so the binary 
invariants for all possible unlabeled graphs span the invariant tensor subspace. 

Theorem 1 The set of binary invariants {B {Q) : Q e G} forms a basis for the DPT 
Hamiltonian coefficient tensors. 
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Proof This result follows from the definition of a basis and that the set of binary invariants 
for all ^ e G is linearly independent and spans the vector space. 

As shown in Reference [20] this result generalizes to any group with any set of tensors 
invariant under that group. 

Therefore, we may resolve the Q tensors at any order as a finite linear combination of 
binary invariants: 

- E ^"°^'(^) [^"''^'(^)]... (18) 

where Guock represents the set of graphs present in the order-0, rank-i? tensor block 
^'rQ^''"'^'' , and the binary invariant B^^°'^^{Q) has the same dimensions as the original Q 
tensor block. The scalar quantity Q'^''°'^^{Q) is the expansion coefficient. 

The resolution of symmetric tensor blocks in the basis of binary invariants in Eq. (18) 
represents a generalization of a technique used at lowest order in Refs. 15 and 21 to arbitrary 
order. This equation also separates the specific interaction dynamics present in Q^^°'^^{Q) 
from the point group symmetry embodied in B^^°^^{Q) . 



IV. THE WAVE FUNCTION THROUGH FIRST ORDER 
A. Lowest-Order Ground State Wavefunction 

Since the Hamiltonian of the lowest-order wave function (Eq. (4)) has the form of a 
A^(A'"-|- l)/2-dimensional coupled harmonic oscillator, the lowest-order wave function will be 
a product of one-dimensional, harmonic-oscillator, normal-mode functions. 

The lowest-order DPT wave function, ^q{c(), for the ground state is given by 

p 

^o(q') = n '^o (v^ ql) , (19) 

where 




There are N{N + l)/2 normal modes and up to N{N + l)/2 distinct frequencies, a number 
which would become impossibly large to solve for if it weren't for the point group 
symmetry expressed in the invariance of the F and G tensors, and the small, A'"-independent 
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number of binary invariants spanning the invariant tensor spaces. In Refs. 15, 21 and 22 we 
have used this Sn point group symmetry to derive both the frequencies and normal modes 
of the lowest-order, Jacobian-weighted wave function for arbitrary N . This analysis results 
in only five distinct frequencies, associated with center-of-mass and breathing modes, radial 
and angular singly-excited state modes, and phonon modes. Each of these frequencies is 
associated with a set of normal modes which transforms under an irreducible representation 
of the Sn point group. 

B. First-Order Wave Function 

Using the Sn point group symmetry expressed in the invariance of the F and G tensors, 
and the small, A^-independent number of binary invariants spanning the invariant tensor 
spaces, in Paper 1 we have also derived the first-order correction to the lowest-order, har- 
monic wavef unction. If we write 

$(q') = (l + 5U)$o(q') + C'(5), (21) 

then A satisfies the commutator eigenvalue equation 

[A,^o]$o = ^i$o. (22) 

To solve this equation, we note that since $o(q') is a Gaussian function, the derivatives in Hi 
and Ho written in normal coordinates "bring down" normal coordinates from the exponent 
so that Hi effectively becomes a 3rd-order polynomial of only odd powers in q'. Then from 
Eq. (22) A is a cubic polynomial and of only odd powers in the normal modes. When A is 
re expressed in terms of internal displacement coordinates, r' and 7' , it is cubic and of only 
odd powers in these internal displacement coordinates. 

The ground-state wave function is also scalar under S^, and so when it is expressed in 
terms of internal displacement coordinates it involves binary invariants which take powers of 
the internal displacement coordinates and couple them together to produce a scalar under 
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V. A TEST OF THE THEORY: THE HARMONICALLY CONFINED, HARMON- 
ICALLY INTERACTING SYSTEM 



The general theory developed in Paper 1 and and Ref. 15, and briefly reviewed in this 
paper, is extensive, and we test it on a non-trivial, interacting, analytic solvable model: the 
harmonically- interacting system of N particles under harmonic confinement: 

+ • (23) 



-if: 



dr 



i<j 



A. The Wavefunction Through First Order 



In Appendix A we solve the harmonically-confined, harmonically-interacting system of 
N particles exactly for the ground-state wave function (see Eq. (A6)), and then from this 

derive the exact analytic perturbation scries for the A^-body wavefunction (weighted by a 
Jacobian) through first order: 
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Repeated indices i , j , ... imply summation from 1 to A?" , while repeated index pairs 



[ij) etc imply the ordered sum 1 < i < j < N . For example, B(^^-^^) 



' N 



. In the above expressions for Ay' and [y']^ fiy/ y' , we are building up 



11 



the invariant polynomials in and 7^^.^^ using the binary invariants as our building blocks. 
The scalar coefficients, A(^) are (derived in the Appendix) 



- ^ (27) 



A(00) = ^ (28) 



(29) 



A(0) = — (30) 

A(^) = ^6 (TV + 1)7^ (31) 

A(e) ^ A{B + CV) (32) 

A(^) = AiB + C£ + J^) (33) 

A(O-) = ^(^ + c(f + (34) 

A{J_) = AiB + C£) (35) 

A(U) = A(B+^-^-g] (36) 



A(0) = A[B+^ + 2gj (37) 

AQ) ^ a(^B + ^^ (38) 

A(~) = AB (39) 

A(a3) = Aeff+^(Aeff-l) (40) 

A(0 0) = ^ (41) 

A(C>-) = roo (42) 

A(o) ^n{I + J) (43) 

A(J)=W(X-^) (44) 

A(~)=7i. (45) 

In the above equations, we have defined 



A=^l + iVA2 (46) 
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(50) 
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This analytic solution through first order is then compared with the DPT wavefunction 

derived from the general formalism of Paper 1 (sec Eqs. (33) and (117)) and implemented 
in Mathematica[23] code. For the case of the general DPT formalism, from Eqs. (19) and 
(20) we find that 

= V^f2q,V, (61) 

[^q']ui,V2 = ^vi,U2^vi , (62) 

and V is the matrix transforming from the internal displacement coordinate vector y' to the 
normal mode coordinate vector q' . The polynomial A of Eq. (21) is similarly transformed 
from a normal coordinate basis to Ay^ of Eq. (24) in the internal coordinate basis. 

In Tables I-IV we compare the binary invariant coefficients, A(^) , from the general 
theory of Papers 1 and Ref. 15 with the above results derived from the full analytical solution 
above for N — 10, 000 particles and two different interparticle interaction strengths, A . One 
value of A features strongly attractive harmonic interparticle interactions, while the other 
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is for weakly-bound system with repulsive interparticle interactions (negative A) for A just 
above the dissociation threshold at A = —1/y/N. 

In both cases, to within round-off-error determined by the machine precision, exact agree- 
ment is found, confirming the correctness of the general formalism of Paper 1, and its im- 
plementation in Mathematica[23] coding. 

VI. SUMMARY AND CONCLUSIONS 

In this paper we performed the first test of a general formalism from Paper 1 for a fully- 
interacting A'"-body wave function through first order in a perturbation expansion. This 
formalism was verified by comparison to a fully-interacting, exactly-solvable model problem. 

The resources required for a solution to a general A^-body problem are understood to scale 
at least exponentially with N , making it very challenging to solve for large- iV systems. [2, 3] 
The present perturbation series will scale exponentially in N if summed to all orders. How- 
ever at first order, the number of terms scale as N^, a scaling, while greatly improved, still 
remains challenging. Nonetheless, this scaling is tamed by expanding the perturbation 
series about a point where the A'^-body system has a highly symmetric structure. In the 
process, the A^-scaling aspect also effectively separates away from the rest of the physics 
allowing the N scaling to be treated as a straight mathematical problem. 

This highly symmetric structure for arbitrary N is obtained as the number of spatial 
dimensions approaches infinity, resulting in a configuration whose point group is isomorphic 
to the (S'jv group. All terms in the perturbation series for the Hamiltonian are then invariant 
under the A^! elements of the group, allowing an expansion in a basis that is also invariant 
under these A^! operations. This restriction results in a comparatively small basis at each 
order which is independent of A^. There arc only seven binary invariants at lowest order for 
any value of N, and twenty-three at next order independent of A^ (except when N is quite 
small when the number is even lower). In this paper we demonstrated the completeness of 
this basis at all orders. Thus order-by-order the wave function, along with other properties, 
may be derived essentially analytically. 

Since the perturbation parameter is the dimensionality of space and not the interaction, 
this approach is equally applicable to weakly-interacting systems for which the mean-field 
approach is valid, and, perhaps more interestingly, strongly-interacting systems for which 
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the mean-field approach breaks down. Paper 1 extends previous work[15, 16, 21, 22] which 
derived energies, frequencies, normal-mode coordinates, wave functions, and density profiles 
at lowest order for quantum systems of confined, interacting particles of any number, . 

The general formalism set forth in Paper 1 for the wave function through first order, 
while essentially analytic, has many moving parts, and so this paper sets out to verify the 
formalism by applying it to the harmonically-confined system of particles interacting via 
harmonic potentials which may be attractive or repulsive. This system is exactly soluble 
in D dimensions, from which we directly derived the dimensional expansion for the wave 
function through first order in terms of the binary invariants. This expansion, directly from 
the exact wave function, has been compared with the wave function through first order from 
the general formalism of Paper 1. Since at each order there arc only a finite number of binary 
invariants to consider, there are only a finite number of coefficients to the binary invariant 
terms that have to be compared. Exact agreement is found between the coefficients obtained 
directly from the analj^ic solution and those derived using the general formalism of Paper 1, 
confirming this general formalism. 

The formalism of Paper 1 is not limited to the harmonic interactions discussed in this 
paper, and indeed is quite general. In previous papers we have examined other potentials, 
in particular the hard-sphere potential in relation to Bose-Einstein condensates, at lowest 
order. While the lowest-order formalism adequately captures the behavior of the system in a 
range of scattering length, a , and particle number N at which the mean field fails, for large 
enough a and/or N the lowest-order wave function no longer has the fiexibility to adequately 
represent the actual system. It is thus desirable to apply the general formalism of Paper 1 
for the wave function, through first-order and verified in this paper, to other systems such 
as the Bose-Einstein condensate for interaction strengths and particle number at which the 
mean field breaks down. 

In principle any observable quantity can be obtained from the wave function and as an 
illustration in Ref. 16 we derived the density profile at lowest order from the lowest-order 
wave function. With the next-order wave function available from the formalism of Paper 1 
we can now derive any observable quantity, such as the density profile, to next order in the 
perturbation theory. 

It is also important to note that while Paper 1 derives the A^-particle wave function to 
next order in perturbation theory, the same basic approach can, in principle, be used to 
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derive yet higher-order terms in the perturbation series. 



APPENDIX A: CONFINED, HARMONICALLY INTERACTING, ANALYTI- 
CALLY SOLVABLE MODEL SYSTEM 

In this appendix wc derive the exact ground-state wave function for a harmonically- 
confined, harmonically-interacting system of particles in D dimensions, and from it derive 
the wave function through first order in 5^/"^ exactly, where 5 — 1/D . 

The Hamiltonian of the harmonically interacting model system of identical particles is 



N 



N 



Making the orthogonal transformation to center of mass and Jacobi coordinates 
where 1 < i < N — 1 , the Hamiltonian becomes 



the sum of N , D-dimensional harmonic-oscillator Hamiltonians, where 

cuint = -y^^^lTA^. (A4) 

Notice two things about the Hamiltonian: it is separable and each component has the 
form of a £)-dimensional harmonic oscillator. Therefore the ground-state solution to the 
wave function in the Schrodinger equation 

H^^E^ (A5) 

is the product of harmonic-oscillator wavefunctions 

N-l 

^{R, {pi}; D) = il^{R- ujt, D) [] V^A; ^int, D) , (A6) 

i=l 

where il^^pi] uj-mt-i D) is the D-dimensional, harmonic-oscillator, ground-state wave function 



r(f) 
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satisfying the normalization condition 

[iljir; uj, D)]'^r^-Ur = 1. (A8) 



The Jacobian-weighted, L — wave function is obtained by folding into the wave- 
function, the square root of that portion of the Jacobian which depends on the internal 
coordinates, i.e. the square root of 



N 



r(D-Ar-i)/2-Q^(D-i)^ (A9) 
where F is the Grammian determinant, so that 

AT iV-l 

^^,J = ^r^D-N-m ^(^-d/^ -q ^.^^^ ^ (aio) 
j=i i=i 

where jV" is a normalization constant ensuring that 

/ {ijih n n ^^^-^ = 1 ■ (All) 

i j<k 



1. A perturbation series in l/y/D for the exact wavefunction 

a. Dimensional scaling 

Now consider transforming to dimensionally-scaled oscillator coordinates 

r ^ D^atr p = D^atp R = D^atR 

_ D^f _ D'^p _ D^R 

\/uh y/^t \f^t 

D^/^f D^l^p D^'^R 



(A12) 



where the cj^ = , and at is the dimensionally-scaled oscillator length of the trap. Both f 
and p are dimensionless. From Eq. (AlO) we obtain 

exp i-^R' 



[(AD)f]V nexp(-:^pl), (A13) 



where 



A = ^ . (A14) 
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b. The large- dimension limit 



To test the general formalism of Paper 1 we need to expand Eq. (A13) about the large- 
dimension limit through first order in S^^^ . In the large-dimension limit the system localizes 
about a structure where all the radii are equal to and angle cosines are equal to 700- To 
derive foo and 700 one applies the condition 



(9^ 



D=oo ^7jfe 
In this endeavor the following results are useful: 



= 0. 



D=oo 



dT 



\D=oo 



D=oo 



\D=oo 



D=oo 



(1 + (A^- 1)700) (l-Too^-^ 
-27oo(l - 7oo)^-' 
fl(l + (N- 1)7^) 



(1 + (TV - 1)700) 



N 



D=oo 
N-l 



2^00 

N 

N 



i=l 



N-l 



D=oo 



D=oo 



D=oo 



(7V-l)f^(l-^^) = (iV-l)po 
2(A^- 1; 



N 



-Tooi^ - 7o 



^Prom Eq. (A15) we obtain the parameters Too and 700 

(A-l)_ 
1)) 

A^+(A-1) 



Too — 

^2 



(7V+(A-1)) 
1 



(A15) 



(A16) 
(A17) 
(A18) 
(A19) 

(A20) 
(A21) 

(A22) 
(A23) 
(A24) 



(A25) 
(A26) 



°° 2(1 + (A^ - l)7oo) 2AA^ 

Equations (A25) and (A26) define the D — > 00 structure about which the system oscillates 
at finite dimension. 
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c. A series expansion about the large-D limit 



To derive the wave function through order S^^^ we perform a series expansion of each of 
the £)-dependent terms in Eq. (A13). 



1 2— exp(f^ 



r (f ) 



N 

n 

i=l 



D-l 



N(D-l) 



exp 2^ 



,i=l 



2f. 



exp 



JV 



:/2 



We also have 



where 



and 



P2 



X 



1 + ^E'^ 



i=l 



/3 ^/ 



3r^ 



3^) +0(5) 



Ea' 

i=l 



R' 


D^oo 


= Rio 


N-l 












i=l 


D^oo 





2f 



oo 



N 

^1/2 / N 



(N N \ 

(l + (iV-l)7oo)E^1+ 
1=1 i<j=l / 

/TV AT AT 



rl/2 / A/ N N 

+ 1^ E(^"^)' + 270O E ^1^1 + E (^1 + ^1)% 

\i=l i<j=l i<j=l J 

N 



+ E 



i<i=l 



+ 



AT 

(51/2 



(A27) 



(A28) 



(A29) 
(A30) 



(A31) 
(A32) 



(A33) 



(N N \ 

(^^-i)(i-7oo)E^i- E 
1=1 i<j=l / 

:l/2 / ^ ^ ^ \ 

— (A^ - 1) J2ir'f - 270O E ^1^1 - E (^1 + ^1)%- 

\ i=l i<j=l i<j=l / 

o A^ 



j<j=l 
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so that 



D 



M-^^'] n 



7V-1 



exp 



— exp 



Pi 



N 



D2f' 



exp 



D2r 



^2 



N 



i=l 



N 



(1 - ^) E % • 



i<j=i 



X exp — ; 



A- 



A- 1 

N 



^ i=l «<j=l 



i<j=l / / 



X 1 + 5^ 



A-1 

N 



N 



i<j=l 



The final bit of the puzzle in the dimensional expansion of Eq. (A13) is the dimensional 
expansion of ]7(-D-Af-i)/4 _ p^j, ^j^jg need Eqs. (A16), (A17), and 



= 



d-fij d-fjk d-fim 



dlij d-fjk d-fki 



d-fij d-fjk d-fji 



dlij d-fjk d-fik 



9ll d-fjk 



D=oo 



D=oo 



D=oo 



D=oo 



D=oo 



D=oo 



D=oo 



D=oo 



D=oo 



D=oo 



27oo(l-7oo)^-' 

-2(l + (iV-3)7oo)(l-7oo)'^-' 





-27oo(l-7oo)^-' 


2(l + (7V-4)7oo)(l-7oor-' 



47oo(l - 7o 



\N-4 



= 



(A36) 
(A37) 
(ASS) 
(A39) 
(A40) 
(A41) 
(A42) 
(A43) 
(A44) 
(A45) 
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0, 



(A46) 



D=oo 



from which we obtain 



^iD-N-m ^ ((1 _ ^^)iV-l (1 + (AT - 1)^^))- 

5i 



X 1 + 



87; 



3 

00 



12(1 -7o,)(l + (iV- 1)7,,) V (l-7oo)=^(l + (iV-l)7oo)' 
6700 



(l-7oo)2(l + (iV-l)7oo) 
1 



[Bi-Wm + 3)7oo)S(0) - 7ooS( J)}7'7' 



+ 



(l-7oo)^ 



1(1 + (TV - 4)7^)i3(/I) - 7oo5(U) + 27oo5(2)}7W 



+ 6(7V + 1)7^[SH7']^ + 0{5)j exp l^-i^^ 



i (A - l)f, 



AT 
i<j=l 



X exp 



2(l-7„o)2(l + (iV-l)7, 



7^ 

/ 00 



(1 + (AT - 1)^^) 



/1 2 



+ 



(A47) 



where the [-B(^)]iyi,i/2,... ^-^e the binary invariants introduced in Paper 1 (briefly reviewed 
in Appendix III) and Q is the graph labehng the binary invariant. The expression 
B{g)X[X2 X3 is shorthand for [B{g)] 

^'i,v2,va [-^i]i^i[-^2 ]j^2[-^3 \v2 where repeated indices 
z/j are summed over, X' is the f' or 7' vector from Eq. (6), hkewise for B{g)X[X2' and 
B{g)X[ . Using Eqs. (A27), (A28), (A35), and (A47), along with 

B{^) ® 5(-) = B{0)+ BU) + Bi^) (A48) 
B{^) ^ B{^) ® B{^) = B{o) + B{^) + B{<:y^)+B{J_)+B{U)) 

+S(0) + S(n) + S(^) (A49) 



B{^)®BiO) = B{0) + ^^^K ^^-^ 



i?(^)( 
in Eq. (A13), with 



3 3 

BUI ^ ^(o-) I B(A) I -5(0^) I B(u) I -^Q) 

2 3 2 2 3 6 



(A50) 
(A51) 



N(D-l) 

foe ' ((l-7oo)^-Ml + (^-l)7oo))^ 



- + 0(5) , 



(A52) 
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we obtain the Jacobian-weighted A^-body wavefunction in Eq. (25) for a system of identical 
particles under harmonic confinement with harmonic interactions. 
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g 


A[A{g)] 


99 


7.0 X 10^1^ 


oo 


2.3 X 10"^^ 


o 


-4.1 X 10-16 


A 


-1.3 X 10-16 


o- 


-6.1 X 10-16 


J- 


1.2 X 10-16 





A [A(g)] 


M 

i — 1 


-1.5 X 10-16 


O 


-6.0 X 10-13 




3.7 X 10-16 




-2.1 X 10-13 





-8.4 X 10-^ 




8.4 X 10-1° 




{AanalyticiG) 



Q 


A[A(g)] 


OO 


-3.6 X 10^16 


00 


-5.1 X 10^11 




9.6 X 10-15 


o 


-3.7 X 10-16 


J 


-2.2 X 10-16 




2.5 X 10-1"^ 



analytic and DPT rank-three, rank-two, and rank-one binary invariant coefficients when N 
10,000 and A = 10. 



Q 


HQ) 


go 


2.3 X 10-11 


000 


1.1 X 10-11 




2.4 X 10-15 




2.2 X 10-15 


oo 


2.0 X 10-15 


qp 


8.2 X 10-16 


o — 


-1.2 X 10-1"^ 



Q 


HQ) 


CO 


-2.9 X 10-15 


oo 


7.8 X 10^19 




4.1 X 10-1^ 


k) 


9.8 X 10-1^ 


LO 


1.0 X 10-19 




6.2 X 10-20 


"ZO 


2.1 X 10-20 



TABLE II: Rank-three and rank-two binary invariant coefficient, AupxiQ) , from the general 
Mathematica code when N = 10, 000 and A = 10. All of these coefficients are exactly zero in the 
exactly soluble analytic solution. 
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A[A{g)] 
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A[A{g)] 


99 


1.6 X 10^1^ 
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-2.0 X 10-8 


oo 


-3.8 X 10"^ 
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1.7 X 10-13 


o 


2.5 X 10"^ 




2.1 X 10-11 


A 


-3.8 X 10"^ 




-3.0 X lO-i'^ 




2.5 X 10-8 





1.7 X 10-3 


J- 


2.2 X 10-5 




-1.7 X 10-10 



g 


A[A{g)] 


oo 


-2.6 X 10^11 


00 


-1.1 X 10-9 




6.8 X 10-13 


o 


-1.6 X 10-'^ 


J 


1.6 X 10-10 




-3.8 X 10-13 



TABLE III: Fractional difference, AA{g) = {AanaiyUc{9) - AopTig)) / AanaiyUc{g) , between the 
analytic and DPT rank-three, rank-two and rank-one binary invariant coefficients when N = 10, 000 
and = -1/10, 000 + 10-i°. 



g 


A(g) 


8^ 


-4.8 X 10-20 


000 


7.9 X 10-2^ 


8^ 


1.9 X 10-12 


o^o 


-1.9 X 10-1^ 


oo 


-2.6 X 10-19 


00 


3.8 X 10-20 


O " 


8.1 X 10-20 



g 


A{g) 


CO 


2.8 X 10-10 


oo 


-1.7 X 10-1^ 


XL 


-2.8 X 10-1^ 


k) 


5.3 X lO-i'' 


LO 


-1.2 X 10-1^ 


-o 


-5.0 X 10-18 




3.3 X 10-20 



TABLE IV: Rank-three and rank-two binary invariant coefficient, A^pxig) , from the general 
Mathematica code when AT = 10,000 and X'^ = -1/10,000 + 10'^^. All of these coefficients 
exactly zero in the exactly soluble analytic solution. 
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